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ABSTRACT 

In light of the growing interest in searching for low mass, rocky planets, we investigate the 
impact of starspots on radial velocity searches for earth-mass planets in orbit about M dwarf 
stars. Since new surveys targeting M dwarfs will likely be carried out at infrared wavelengths, 
a comparison between V and Y band starspot induced jitter is made, indicating a reduction 
of up to an order of magnitude when observing in the Y band. The exact reduction in jitter 
is dependent on the photosphere to spot contrast ratio, with greater improvements at smaller 
contrasts. 

We extrapolate a model used to describe solar spot distributions to simulate the spot pat- 
terns that we expect to find on M dwarfs. Under the assumption that M dwarfs are near or 
fully convective, we randomly place starspots on the stellar surface, simulating different lev- 
els of spot coverage. Line profiles, distorted by spots are derived and are used to investigate 
the starspot induced jitter. By making assumptions about the degree of spot activity, detection 
limits for earth-mass planets in habitable zones are simulated for between 10 and 500 obser- 
vation epochs. We find that ^ 50 epochs are required to detect 1-2 M 9 planets (with < 1 per 
cent false alarm probability) orbiting slowly rotating 0. 1 and 0.2 M Q stars. This sensitivity de- 
creases when typical rotation velocities and activity levels for each stellar mass/spectral type 
are considered. No detections of below 20 M ffl planets are expected for ^ 500 observations 
for the most active stars with v sin i 20 km s _1 and dark spots. 

Key words: (stars:) planetary systems stars: activity stars: atmospheres stars: spots tech- 
niques: radial velocities 
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1 INTRODUCTION 

Nearly 500 extrasolar planets in over 400 planetary systems have 
been discovered and yet < 6 per cent of those with known or es- 
timated host star masses are M dwarf planetary systems. Since M 
dwarf masses may be as low as one tenth of a solar mass, it is not 
unreasonable to expect that their most mass ive planets might also 
be correspondingly smaller flda & Lin 2005) and thus more diffi- 
cult to detect than the ubiquitous gas giants found orbiting F, G & 
K dwarf stars. On the other hand, planets of a fixed mass orbit- 
ing lower mass stars become easier to detect, owing to the reduced 
mass ratio of the system. Detecting l ow-mass planet s (which are 
predicted in significant numbers by llda & Linl Q2Q05)) is however 
no easy task, since the planet-induced stellar velocity amplitudes 
are small, even for M dwarfs. Additionally, in the case of M dwarfs, 
there is an observational bias against detecting planets at optical 
wavelengths (where most searches have so far been carried out) ow- 
ing to the relatively low fluxes compared with redder wavelengths. 
Observations carried out at near infrared wavelengths, where the 
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host M stars are several magnitudes brighter than at shorter wave- 
lengths, open up the possibility of carrying out larger surveys of 
lower mass stars that are capable of detecting low-mass planets. 

Radial velocity surveys are still by far the most success- 
ful means of detecting planets. Despite the difficulties of ob- 
serving in the infrared, there are a number of fledgling projects 
aim ed at carrying out plane t searches. Both lRamsev et all ( 2008) 
and Steinmetz et 

aD [2008) showed that high precision (< 10 
m s _1 ) radial velocity measurements could be made at short in- 
frared wavelengths via observations of the Sun. Seifahrt & Kaufl 
d2008l) demonstrated that similar precision was possible on short 
timescales using the Cryogenic high-resolution infrared echelle 
spectrograph (CRIRES), operating at L-band wavelengths, at the 
Very Large Telescope (VLT). More recently, an ammo nia gas cell 
has be en employed to achieve ~ 5 m s _1 precision dBean et aU 
2009b) over six month timescal es with CRIRES operating at K- 
band wavelengths. This survey dBean et al.l l2009a) has not found 
evidence for the m assive planet claimed to be orbiting the low-mass 
M dw arf VB 10 dPravdo & Shaklanll2009l ; IZapatero Osorio et all 
l2009h . 

With infrared multi-order spectrographs working at high res- 
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olution (R > 50, 000), precision radial velocities of order 1 m s 
are expected to be achieved in the next few years. Instrumental 
precision of 1 m s" 1 is a desirable goal when we consider that a 
1 M9 planet orbiting in the habitable zo ne of 0.1 Mr star in duces 
a ~ 1.7ms -1 radial velocity amplitude. iRivera et al] d2005h have 
nevertheless already estimated < 1.5 m s" 1 radial velocity jitter 
for the M4 dwarf, GJ 876 (at optical wavelengths). 

A major cause of radial velocity jitter is the presence of mag- 
netic activity induced s tarspots that distort absorption line profiles 
dSaar&Donahuell 19971) . Any time variant asymmetries in line pro- 
files due to starspots will lead to biased measurements of dynami- 
cally induced radial velocities from orbiting bodies such as planets. 
If the radial velocity signal of a planet is of similar magnitude or 
smaller than spot induced jitter, any candidate planet signal may not 
be recovered with reasonable allocations of observing time. In this 
paper, we simulate semi-realistic starspot patterns for M dwarfs. 
By making use of starspot size distribution models for the Sun and 
by considering indirectly observed starspot distribution patterns for 
M dwarf stars, we model the subsequent line profiles from which 
radial velocity jitter is determined via cr oss-correlation. Thi s ap- 
proach extends simi lar work carried out bv lDesort et alj d2007b and 
iReiners et al.l d2010h to more realistic spot distributions that are so- 
lar spot distributions extrapolated to active stars. The met hod is 
more analogous to the recent work by lLagrange et al.l J2010t> that 
investigates the effect that solar sunspot activity would have on de- 
tection thresholds of an earth-like planet orbiting in the habitable 
zone. Here, we use our models to determine the detection thresh- 
olds for earth mass habitable zone planets orbiting M dwarf stars 
that exhibit different equatorial rotation velocities and starspot ac- 
tivity levels. In fj2]we discuss evidence for M dwarf starspot dis- 
tributions and introduce our extrapolated solar model. In |J3] we 
briefly describe the methods used to generate line profiles. A com- 
parison of the relative spot amplitudes for different photosphere- 
to-spot temperature contrast ratios (T p /T 3 ) is made in $3.2\ The 
radial velocities induced by stars that exhibit the spot distributions 
introduced in $\2.2\ are then investigated in «j3 .31 Here, we inves- 
tigate the effects of spots for three stellar masses at extremes of 
the expected T p /T s values. We then carry out detection threshold 
simulations for Earth-mass planets in ij4]before finally considering 
possibilities that generalise our specific model cases further in ij5] 



2 M DWARF SPOT MODELLING 
2.1 Observed starspot distributions 

Compared with earlier spectral types, our knowledge of starspot 
distribution patterns on M dwarfs is less complete. For earlier spec- 
tral types, publications that predict and report observations of these 
distributions are numerous. The solar analogue is our chief refer- 
ence point for which we observe spots appearing chiefly at 0° - 40° 
latitudes. There is also considerable evidence that other stars ex- 
hibit similar starspot ac tivity to the Sun. Obs ervations made by the 
Mount Wilson Survey dBaliunas et alJI 19951) . which has observed 
changes in the rotatio n periods of a number of stars over decades 
dDonahue et al.| [l996). suggests that similar magnetic dynamo pro- 
cesses are at work. The appearance of spots within defined (low) 
latitude bands has been attributed to the interface dynamo process 
by which the magnetic fields responsible for the spots are gener- 
ated. A magnetic dynamo located at the boundary between the ra- 
diative core and convection zone together with radial transport of 
flux through convection, is able to expl ain the appearance of photo - 
spheric flux at low-mid latitudes only ( Moreno-I nsertis et al J 1992b . 



For more rapidly rotating stars which can be in directly imaged 
via th e Doppler imaging process, this scenario dSchiissler et al.l 
1 1996b is less succ essful at predicting the starspot distributions 
l Barnes et al Jl 1998b . However spots still often appear within spe- 
cific latitude regions in Doppler images. 

Unlike G and K stars, owing largely to their inherent faint- 
ness, M dwarfs are less well studied. By mid-M, it is predicted 
that stars become fully convective and thus the standard aQ, inter- 
face dynamo process can no longer operate. There is nevertheless 
a great deal of evidence indicating that fu lly convective M dwarf 
stars are magnetically active. A study bv lGiziset all f 2000) found 
that the percentage of obje cts with Ha appearing in emission peaks 
at late spectral type while Mohantv & BasrJ d2003l) found that the 
equatorial rotation velocity at which Ha activity saturates appears 
at higher values for mid-late M spectral types than for early-mid M 
spectr al types. This was however later ruled out bv lReiners & Basril 
(2010). The topic is too large to review at length here where we are 
concerned primarily with starspot distributions. 

Observations of rapidly ro tating early M-dwarfs via indi- 
rect Doppler imag ing techniques ( Barnes & Collier Cameronl200ll : 
iBarnes et alj|2004l) reveal that they are more uniformly co vered 
with spots t han e arlier spectral types. iDonati e t al. ( 2008]) and 
Mori n et al. I J2008b have found, via magnetic Doppler imaging 
using Stokes V measurements, that the magnetic topology of M 
dwarfs change at approximately spectral type M4V. Although a 
switch to a fully convective dynamo might b e expected to l ead to 
a m ore uniformly spot ted star, the results of IDonati et al.l (2008) 
and lMorin et alj (2008) appear to counter-intuitively indicate that 
fields beco me more dipolar . This p henomenon has also been inves- 
tigated bv lReiners & Basril (2009) who find that in fact more than 
85 per cent of the magnetic flux is stored in magnetic fields that 
are invisible to Stokes V. However, using all Stokes components is 
necessary for a full description of a star's magnetic topology. Ad- 
ditionally, the field from the darkest magnetic regions (e.g. spots) 
is not visible in Stokes V owing to the large contrast seen at opti- 
cal wav elengths. The high average fields of a few kG found on M 
dwarfs djohns-Krull & Valentilll996l : lReiners & Basrill2007l . l2010b 
are unlikely to be concentrated in small spots because very high lo- 
cal fields would be necessary. They are perhaps evenly distributed 
across the surfaces of the stars. The traditional picture of spots as 
the only concentrations of magnetic flux may therefore be mislead- 
ing in these stars. It is thus clear that neither magnetic imaging nor 
surface brightness imaging at optical wavelengths are able to give 
a complete picture of the starspot/magnetic field topology of stars. 

While starspot patterns (e.g. the latitudes at which spots 
appear) on more rapidly rotating G and K stars may vary as 
a function of rotation, it is not clear whether such changes 
take place among fully convective M dwarf stars that are ex- 
pected to generate magnetic fields via a turbulent dynamo pro- 
cess. Hence those surface brightness images derived for M dwarfs 
dBarnes & Collier Cameronll200ll : IBarnes et al.ll2004b may also be 
representative of slower rotators. Moreover, at la ter spectral types, 
M7-M9.5 for example, IReiners & Basril d2010b find no correla- 
tion between rotation and magnetic flux generation. Further ev- 
idence that starspots may be more uniformly distributed comes 
fro m observations of red uced starspot induced lightcurve variabil- 
ity. iMessina et alj d2003b have shown that already by spectral type 
K6 - M4, the maximum starspot induced photometric variability is 
around a factor of 2 lower. Other individual studies o f mid-M and 
late-M (M5 & M9) dwarfs dRockenfeller et alj|2 006a b) also show 
peak-to-peak amplitudes of order 0.05 or less in G, R and I pho- 
tometric bands. A more uniform distribution of spots (rather than 
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Figure 1. The distribution of spots for Models 1 (top left) to 6 (bottom right) for an M dwarf. Models 1 and 2 are analogous to solar min and solar max activity 
levels while Model 3 represents a high solar activity case. Models 4 to 6 are included for completeness (see Table 1). 



clustering of spots in one or two active regions) would be one ex- 
planation for the observed reduction in photometric induced light 
variations in M dwarf stars. 

A series of investigations into spot coverage factors 
jNeff et alJ[l995l ; l0'Neal et al Jl 19961. 1998 & 2004) have been car- 
ried out using TiO as a tracer of cooler temperatures on G & K 
stars of different spectral class, including dwarfs. These studies, 
indicate typical spot coverage of 20 - 50 per cent for active stars. 
This may seem surprising, but indicates that lightcurve analyses 
and Doppler imaging studies (that typically find of order 10 per 
cent spot coverage) are not sensitive to an underlying spot distribu- 
tion. It is again unclear whether these spots arise from a boundary 
dynamo or whether a turbulent dynamo is responsible. Based on 
these findings, investigations of high levels of spot coverage seem 
warranted. A threshold level of spot coverage may be expected (for 
uniform spot coverage) at which any radial velocity jitter effects no 
longer increase in magnitude. 

With a largely unknown starspot pattern in moderately rotat- 
ing M dwarfs and the general prediction that a distributed dynamo 
should produce randomly distributed spots, we carry out a num- 
ber of simulations to assess the detectability of planets around M 
stars. Without strong evidence to the contrary, we assume that more 
rapidly rotating stars are more spotted. Synthetic starspot models 
are used to generate line profiles and to investigate the radial veloc- 
ity amplitudes resulting from non-uniform line profiles. These ra- 
dial velocities are then used to determine our ability to detect low- 
mass planets which are in habitable zone orbits around M dwarf 
stars. 



2.2 Modelling randomly distributed spots 

We synthesize spot maps which follow a log-normal size distri- 
bution on the surface of an immaculate star. T he Doppler imag- 
ing co de "Doppler Tomography of Stars" (DoTS) dCollier Cameron! 
1 1997b was then used to place the spots on the surface of a model 
star. The input parameters to DoTS for modelling spots were first 



introduced in |jeffersU2005h . They are, where x is a random number 
in the range (0 ^ x ^ 1): 

(i) longitude: randomly distributed between 0° and 360° 

(ii) latitude: — ~ < 6 < |, following 8 = arcsin (2x + 1) with 
O^x^l, to eliminate an artificial concentration of spots at the pole 

(iii) spot radius: computed using the previously described log- 
normal distribution as tabulated in Table 1 

(iv) spot brightness & spot sh arpness: mod elled to obtain an 
umbral to penumbral ratio of 1:3 (Solanki 1999) 

The key difference when compared with the solar case is that 
spots are allowed at all latitudes. The sola r spot size distribution has 
been determined by Bogd arTet alj dl988l) from direct observations 
taken from the Mount Wilson white- l ight p late collection covering 
the period 1917-1982. iBogdan et ail dl988l) show that the number 
of sunspots, N, with umbral area, A, is given by 



dN 



dN\ 
dA J 



exp 



(In A -In (A)) 2 
2 In a a 



(1) 



where the constants (A) and a a are the mean and geometric stan- 
dard deviation of the log-normal distribution of the sunspot areas, 
and (4gr) is the maximum value reached by the distribution. 

V dA I max J 

For the case of the Sun, these values are tabulated in Table 1 where 
model 1 is for an inactive Sun and model 2 is for an active Sun. 
Cool, young potentially planet hosting stars are expected to be more 
active than the Sun. For these stars we use the extrapolations of 
ISolankil Jl999h als o show n in Table 1 . All starspots are modelled, 
following Solank| dl999h with circular umbral areas. A starspot is 
defined such that the umbral region is at the temperatur e of the 
spot (t abulated in Table 2). Following the assumptions of Solanki 
(1999), we also include penumbral regions with umbral to penum- 
bral areas of 1 :3 (i.e. radii ratios of 1 :2). The penumbral regions are 
defined such that their intensity is equal to half the difference of the 
photosphere and spot intensities. 

We have included the solar values for R' H K at solar min and 
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Model 


1 


2 


3 4 


5 


6 


<T A (xlO- 6 A 1/2Q ) 


3.8 


5.0 


6.8 9.2 


12.2 


15.8 


(dN/dA) maa; 


5 


25 


65 125 


205 


305 


Umbral Spot Coverage 


0.03% 


0.3% 


1.6% 6.1% 


18% 


48% 


Total Spot filling 


0.03% 


0.3% 


1.9% 9.0% 


29.5% 


62.4% 


Solar: 


min 


max 










-5.0 


-4.85 




-> TiO obsns. 



Table 1. Tabulation of the input parameters to the log-no rmal size distribution of star spots given by equation[TJ The parameters derived by Bogdan et al. 1 1988) 
for the Sun are data sets 1 & 2, and those calculated bv lSolanld fl99l for active stars are data sets 3-6 (Ai/ 2 Q = 27ri?g). The spot filling is determined 
directly from DoTS and differs from the umbral spot coverage since we also include penumbral regions. The spot coverage and filling are the s ame for models 
1 & 2 since the spot sizes are close to the pixel resolution in these cases. R'„ y values correspond ing to the solar min and max cases are given iLagrange et alj 
|2010|) and the upper limit for starspot coverage derived from TiO studies (e.g lO'Neal et aljl9 98) is indicated. 



solar max, as derived bylLagrange et al] (feoioh who used t he rela- 
tionships o f lNovesetal.Ul984h and Lockwood et "al] b007h . These 
values also give an indication of the kinds of a ctivity levels that ar e 
routinely found in optical planet survey targets dTinnev et al.l2002h . 
We also indicate the models (i.e. numbers 5-6) that correspond to 
the upper limit of ~ 50 per cent spot filling fraction determi ned in 
active G and K stars fo'Neal et alll998r . lO'Neal et al.l20o4) . Con- 
sideration of such heavily spotted stars is of interest, especially if 
the contrast ratio, T p /T 3 is also much lower in mid-late M stars 
(see g27T]c& ij33](. This may enable more active stars to be included 
in radial velocity surveys. 



3 GENERATING SYNTHETIC RADIAL VELOCITY 
DATA 

The imaging code, DoTS, enables us to generate line profiles us- 
ing a 3D stellar model with arbitrarily placed spots. Radial velocity 
measurements are then made directly from the line profiles. DoTS 
works with two intensity profiles for each image pixel, representing 
the photospheric and spot temperatures. The degree of spot filling 
for each pixel is represented by a value in the range 0.0 - 1 .0, where 

1.0 represents complete filling (i.e. an image pixel at the spot tem- 
perature). 

Relative fluxes for a star with a given temperature are de- 
termined from the abso lute magnitude and r adii determined by 
Bara ffe et al] i 19981) and lChabrieretalJ d2000h for low mass stars 
and substellar objects. Hence we are able to determine the fluxes 
for a star of any mass with a specified model photospheric tem- 
perature and an attributed spot temperature. We use two specific 
photosphere/spot temperature contrast extremes in later sections of 
the paper. However, in the following section we begin by examin- 
ing different photosphere/spot contrast ratios. Thr oughout this pa- 
per, we use non-linear limb-darkening values from lClaretl ( feOOOh to 
represent the radial variation of intensity for both photospheric and 
spot temperatures. 

3.1 Choice of near infrared passband 

Our choice of infrared passband for the simulations was determined 
by the observation that some near infrared wavelength regions con- 
tain a greater degree of spectral line information, in the sense that 
more sharp features are found in these regions. This was noted for 
example by Rei ners et alj d2010h who found that the highest radial 
velocity precision was achieved in the Y band when compared with 
the J and H bands. The Y band outperformed the J and H bands in 



terms of achievable spectral precision for all M dwarf spectral types 
(see their Fig. 5). Only for the early to mid M spectral types, owing 
to the much higher spectral information content, does the V band 
still outperform the infrared bands despite much lower flux. As we 
shall see in the next section, the main advantage of observing at 
infrared wavelengths is the reduction in starspot induced jitter. For 
the following simulations, we therefore focus on the Y band which 
appears to offer the best chances of detecting low amplitude signa- 
tures that may arise from the reflex motion of orbiting planets. 

3.2 Relative visible and near infrared spot induced 
amplitudes 

We first investigate the effect of contrast variations on the ampli- 
tude of the radial velocities. Observations made at two different 
wavelengths, the V-band and Y-band regions (centered at 5450 A 
and 10350 A respectively), are simulated. We begin by placing a 
single spot with radius 10° on the equator of a star inclined at 90°. 
The stellar line profile is calculated at a number of phases for which 
the starspot is visible. The only differences between phases are the 
intensity (due to foreshortening angle) and radial velocity of the 
spot contribution to the integrated line profile. Since we are inter- 
ested in the relative amplitude of starspot jitter only in this section, 
the simulated v sin i is not important (we used v sin i = 10km s- 1 ), 
providing it is greater than the instrumental resolution. 

Using the above parameters, we show in Fig.|2](top), the am- 
plitude induced as a function of wavelength for the case where T p = 
3250 K. The curves represent high contrast (T s = 2000 K) and low 
contrast (T s = 3000 K) scenarios. The decrease in amplitude as a 
function of increasing wavelength is more pronounced for the low 
contrast case. Here, a secondary effect, namely the relative equiv- 
alent width of the lines at the photospheric and spot temperatures, 
also becomes important and is responsible for the increase in am- 
plitude at A — 10350 A and A = 13000 A when compared with 
shorter and longer wavelengths. This effect is discussed further be- 
low. The results are in broad agreement with those presented in 
lReinersetaH ( l2010l) . 

In Fig.|2](bottom) we show the relative V-band/Y-band spot in- 
duced amplitudes, Ky /Ky , for various spot contrast ratios. There 
are a number of factors that characterise the magnitude of the 
starspot radial velocity induced jitter in a given passband, or spe- 
cific wavelength range, the most important being (1) the contrast ra- 
tio between starspot and photosphere and (2) the normalised equiv- 
alent width of the local intensity profile. Additional, secondary ef- 
fects come from changes in limb darkening coefficient which we 
also model using non linear limb darkening values (see Fig. 3 of 
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Figure 2. Top: Radial velocity amplitude as a function of wavelength for 
an equatorial spot of radius 10° . The curves are shown for high contrast and 
low contrast scenarios with T v = 3250 K. Bottom: The radial velocity am- 
plitude induced by spots in the visible (V band) compared with the infrared 
(Y band). For the labelled photospheric temperatures in the range 2000 K 
< T p < 4000 K, spot amplitudes ratios are calculated for 1500 K ^ T 3 ^ 



Claret (2000)). The local intensity profile may be considered as the 
local un-broadened spectrum of a non-rotating star. The combina- 
tion of the Doppler shifted local intensity profiles from all points on 
the star, or in our case, all pixels on the model star, are summed up 
to give the observed profile (it is this summation of local intensity 
profiles that is also employed to derive indirectly resolved Doppler 
image s dVogt & Penrod 1 19831 : FVogt et al.l [19871 ; ICollier Cameronl 
l200ll) of spotted stars). The effects of (1) an d (2) are investigated 
and discussed in detail in lReiners et al.l (2010). Here, we summarise 
those effects in order to illustrate the trends seen in Fig. 2. The 
Doppler shifted local intensity corresponding to the location of a 
cool spot on a star is much less compared with the equivalent photo- 
spheric contribution. The corresponding photosphere/spot intensity 
ratio is typically greater than an order of magnitude for standard 
solar-like spots that exhibit T p — T s ~ 1000 K. In this instance, 
because the spot contribution to the Doppler broadened profile is 
so much smaller than the photospheric contribution, the spot pro- 
file EW and shape are relatively insignificant. In other words, spots 
contribute little flux, both in the line and the continuum, when the 
contrast between photosphere and spots is large. Hence the missing 
light at a given Doppler shifted position is relatively independent 
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Figure 3. The Y band rms radial velocity variation for each of the 6 simu- 
lated spot models. Pairs of values (connected by a line) are shown for each 
stellar mass (0. 1 , 0.2 and 0.5 Mq cases are labelled for model 2) and stellar 
v sin i values of 1 , 2, 5 , 1 0, 20 & 50 km s - 1 . Each pair of values represents 
the upper and lower variation for the two extreme F P /F B contrasts listed in 
Table 2. 



of wavelength. This is illustrated in Fig. 2 (bottom) by the obser- 
vation that for any T p , Ky/Ky is close to unity for the smallest 
values of T s . The value of Ky / Ky increases as the spot tempera- 
ture (and hence intensity) increases (i.e. as T p /T s decreases). This 
trend arises simply because the contrast ratio for fixed T p /T s is 
much lower in the infrared than in the optical. The wavelength de- 
pendence of the RV signature is illustrated for blackbody fluxes in 
IReiners et alj bold ) . Here we use fluxes derived fro m the mod- 
els of lBaraffe etall ( fl998l) and Ichabrier et alj (2000) for 5 G yr 
stars. The top-right point on each curve in Fig. 2 (bottom) indi- 
cates T p — T s = 250 K. This small temperature difference leads to 
starspot induced jitter that is several times smaller in the Y band as 
compared with the V band. 

For the T p = 3000 K and T p = 4000 K cases, reductions 
in starspot induced jitter of an order of magnitude are achieved. 
These augmented reductions are due largely to the behaviour of 
the relative line intensities in the spot and photosphere at V and 

Y band wavelengths. While the continuum contrast is low, the rela- 
tive equivalent width and depth of the lines is greater in the V band, 
compared with the Y band, leading to the large observed ratio s in 
Fig. 2. Again, this effect is discussed by IReiners et alj j2O10h in 
§4.1.2 of their paper. In essence, for low contrast ratios, the spot in- 
duced jitter is much smaller since there is little difference between 
the intensity of the photospheric and spot contributions. However, 
the ratio of line depths at T p and T s then becomes important. If 
the continuum contrast ratio (1) is small but the line depth of the 
spectral lines in the spot is greater than in the photosphere, the 
two effects may cancel. In other words, the emission bump created 
by the cooler spot spectr um is cancelled by additional absorption 
due to deeper lines (see IReiners et al. I < |2010|) . Fig. 11). We de- 
ri ve a mean line depth for ea ch passband region from the models 
of lBrott&Hauschildtl d2005h - For T p = 3250 K, we see in Fig. 2 
(bottom) that Ky/Ky is much smaller at low T p /T 3 . This arises 
because of a significant growth in equivalent width of lines in the 

Y band below this temperature. For the cases where the continuum 
contrast ratio is small, the centre of the lines are then much deeper, 
leading to higher contrast in the lines. As we shall see in ij3.3l this 
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leads to larger starspot induced velocity amplitudes for our simu- 
lated 0.2 Mq star in low contrast scenarios. 

Our results in Fig . 2 are in close agreement with those pre- 
sented in Reine rs et al.l fcOlOl) . A marked difference between the 
change in radial velocity induced jitter as a function of wavelength 
is seen when moving from their "toy" model results of Fig. 1 1 to 
the model atmosphere result s (analogous to our re sults) of Fig. 12. 
When comparing K v /K Y , IReiners et all feo id) find an approx- 
imate decrease in the starspot induced jitter of ~2 times & ~5 
times for T p = 3700 K and T p = 2800 K respectively, with 
Tp — T s = 200 K. We similarly find starspot induced jitter ratios of 
~ 5 & ~ 6 for T p = 3750 K and T p = 2750 K respectively, with 
Tp — Ts = 250 K. The discrepancies in the results likely arise from 
differences in the models, the line depths used to represent the spot 
and photosphere local intensities and to a smaller degree the dif- 
ferences in limb darkening models used. Nevertheless, the results 
are in broadly close agreement, indicating a reduction in radial ve- 
locity induced jitter in the Y band over the V band of order 5 for 
T p - T 3 = 250 K. 

Fig. [2] shows a clear advantage of observing at near infrared 
wavelengths, when the contrast ratio between spots and photo- 
sphere is not large. At higher contrasts, the advantage is less ob- 
vious. Obtaining an estimate for reasonable values of T s vs T p is 
not straightforward, but clearly has important implications for fur- 
ther simulations. This issue is considered further in the following 
section. 

3.3 Radial velocities induced by random spot distributions 

The absolute radial velocity amplitudes depend not only on the 
starspot distri butions but also on th e photosphere/spot contrast ra- 
tio. Following IReiners et al.l d2010h . we have investigated two ex- 
treme scenarios for M dwarf stars with three different masses. In 
Table 2, we list the chief parameters for the simulated stars, includ- 
ing the photospheric and spot temperatures. 

Stars with axial inclination, i = 90° are simulated according 
to the parameters tabulated in Table 2 and for each of the 6 stellar 
models described above. A total of 36 line profiles (i.e. every 10°) 
with a range of v sin i values of 1, 2, 5, 10, 20 & 50 km s" 1 are 
generated for a complete stellar rotation period. The local intensity 
profile used to represent the rotation profile of a non-rotating star 
is a synthetic Voigt profile with an instrumental resolution that is 
appropriate (R = 70, 000) for high resolution IR spectroscopy (e.g. 
Ijones et alj|20091) . 

Figure [3] plots the starspot induced rms radial velocity varia- 
tion for each spot model (introduced in §2.2) and v sin i value. The 
pairs of values for each stellar mass of 0.1, 0.2 and 0.5 (points con- 
nected by lines in Fig. 3) represent the T S1 = 0.65T P (greater rms 
RV) and T 32 = T p - 200 K (lesser rms RV). As expected, there is 
little variation due to the stellar mass for the T S1 — 0.65T P point 
since the relative contrast ratio is unchanged. The lower point in 
each case shows more variation for a given model since the ratio 
T S2 = T p - 200 K is not fixed. Moreover, for the T p /T s combi- 
nation of 3250 K / 3050 K, we see that the radial velocity variation 
is between 1.5 and 6 times larger than for the corresponding 3750 
K / 3550 K case. This difference is due to the effects discussed in 
the §3.2, where the equivalent width of the lines grows for tem- 
peratures close to, but below 3250 K. The large range in factors 
of between 1.5-6 is likely due to the spot sizes and number that 
define each model and their relationship with the Doppler/spatial 
resolution limit which improves (i.e. smaller spots resolved) with 
increasing v sin i. This will ensure that for low contrast spots on 



our 0.2 Mq model, the starspot jitter will be greater than for the 
0.1 Mq and 0.5 M models. 

For model 1, with vsini = 2 kms (i.e. a slow rotator anal- 
ogous to the Sun) the approximate range of velocities (consid- 
ering all three stellar masses) for the two T p /T s extremes are 
~2-6 cms" 1 , while for model 2 the range of observed values is 
~ 4 cms -1 - 27 ms" 1 . This is an or der of magnitude smal ler than 
the values observed for the Sun by lLagrange et "all feOld) which 
show equivalent radial velocity jitter of solar minimum and max- 
imum periods of up to 0.6 m s _1 and 2 m s _1 respectively. The 
reduction in contrast ratio in the optical vs the infrared, and the fact 
that spots may appear at all latitudes, act to reduce the simulated 
radial velocity jitter for models 1 and 2. This latter point is impor- 
tant since placing only a few small spots randomly at all latitudes, 
as opposed to restricting them to solar-like low-latitude bands may 
lead to significant reductions in the measured jitter (i.e. spots at 
low latitude induce the greatest jitter). By the time a star reaches 
vsini = 10 km s" 1 , we find rms velocities of ~ 4-38 cms" 1 and 
~ 17 cms -1 - 1.9 m s" 1 for models 1 and 2 respectively. Clearly, 
the solar minimum and solar maximum analogue cases (models 1 
and 2) will enable the greatest possibility of detecting lower-mass 
planets since they yield the lowest rms RV variations. As expected, 
the rms jitter peaks at Model 6 which corresponds to a star with a 
62.4 per cent starspot filling fraction. 

For those early M dwarfs t hat are fast rotators, the spot sizes 
derived from Do ppler images (Barn es~& Collier CameroiJ 1200 it 
iBarnes et alJl2004D are much larger than is seen for the Sun. It is 
not clear however whether the spots seen in Doppler images are 
actual individual spots or indeed unresolved spot groups. Indeed if 
we take Model 6 and reconstruct the surface spot distribution for a 
moderately rotating case, we derive images r e sembling the Dopple r 
images of lBarnes & Collier Cameror] d200ll) : lBarnes et all d2004l) . 
We have also investigated the effect of larger spots on the radial 
velocities by scaling the sizes used for models 1 - 6 by factors of 2, 

5 & 10. In other words, a model identical to each of those shown 
in Fig. [T] was created, but with all spots scaled up in radius. The 
starspot induced rms jitter maximum in Fig.[3]is found to increase 
by a factor of ~ 2 for a spot size scale factor of 10. The main effect 
however is a change in the model number at which the peak jitter 
occurs. For scale factor 2, the rms jitter is less peaked for model 

6 when compared with model 5. For scale factor 5, a peak occurs 
between models 3 & 4, while the peak is at model 3 by scale factor 
10. However, even for the more modest scale factor 2 scenario, the 
planet detection thresholds (see §4) for model 2 will more closely 
resemble those of model 3 where the spot scale factor was unal- 
tered. In the following sections, we do not scale the spot sizes, but 
use the spot size distributions as previously defined in §2.2 

3.4 Stellar mass vs v sin i and spot filling 

The effect of v sin i is an important factor that must be considered 
when determining the RV jitter. While Fig. 3 shows that model 1 
does not yield jitter of > 1 m s" 1 , even at v sini = 20 km s - , a 
very active star such as that represented by Model 6 only yields jit- 
ter of order 1 m s _1 or less for the low contrast case when iisini 
= 1 km s _ . There is approximately a 1:1 correspondence between 
jitter and v sin i for models 3 - 6. In other words, doubling the v sin i 
doubles the rms jitter. However, for models 1-2, the increment is 
small owing to the finite spectral resolution of 70,000 which corre- 
sponds to 4.3 km s _1 . 

Since the rotation-activity relation seen at earlier spectral 
types is also seen to persist in M dwarfs, through observation of Ha 
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Table 2. Simulated stellar mass, photospheric temperature, T v , spot temperatures T S1 = 0.65T P (to the nearest 50 K) and T S2 = T p — 200 K. Also tabulated 
are the periods of orbiting habitable zone planets for each stellar mass and the stellar radial velocity amplitudes induced by orbiting planets of mass 1, 2, 5, 10 
&20 M e . 




Figure 4. Example of a simulated stellar RV curve for a 2 Me planet orbit- 
ing a 0.2 Mq star (red/solid line). The RV curve is shown with added jitter 
for an instrumental resolution (IP) of 2 km s —1 , stellar v sini = 5 km s" 1 , 
and starspot model 2. 

( Moha ntv & Basrill2003h . it seems reasonable to assume that the 
slower rotators will be the least active stars while the more rapid ro- 
tators exhibit higher degrees of spot coverage. Unfortunately, since 
there have been no successful studies of spottedness of M dwarfs 
dO'Neal et alj|2005l) . relating «sini to the spot filling fraction is 
not possible. We must therefore make assumptions about the de- 
gree of spottedness as a function of rotation velocity. Based on so- 
lar ob servations and on the spot filling fraction measured on othe r 
stars dO'Neal et alJl 19961 ; lO'Neal et al.lll998l ; b'Neal et al]|2004h . 
we assume that models 1-2 are appropriate for v sini = 2 km s _1 , 
and models 3-6 are appro priate for v sin i = 5 - 50 km s~ . 

I Jenkins et alj d2009h have recently determined the v sin i val- 
ues for a number of M dwarfs, thereby extending the sample 
of measured rotation values. This study reinforced the observa- 
tion |Delfosse et al.ll99 8; Mohanty & Basri l2003l ; lReiners & Basril 
2008) that amongst field M dwarfs, the earlier spectral types ex- 
hibit slower rotation compared with mid and late M spectral types. 
In other words, while earlier spectral types are more likely to be 
slower rotators, enabling more precise radial velocities to be deter- 
mined, the higher mass of the star sets a lower limit on the reflex 
motion due to an orbiting planet of given mass. Conversely, the de- 
tectable mass limit for an orbiting planet is adversely affected by 
the more moderate rotation of mid to late M dwarf stars, while the 
lower mass is more favourable for detecting low mass planets. 

For sp ectral types M2V an d M4V (0.5 M Q and 0.2 M re- 
spectively), Jen kins et alj J2009I) (see their Fig. 9) find the median 
ii sin i is 3-5 km s _1 whereas for spectral types M6V (0.1 Mq), 



vsini ~8 km s" 1 . However, velocities from 1 km s -1 to several 
10s of km s _1 are found for all stars in the spectral range consid- 
ered, albeit with fewer slow rotators found at later spectral types 
(i.e. for masses ^0.1 Mq (M6V)). We therefore simulate a range 
of rotation velocities for each spectral type, but urge the reader to 
bear in mind the trend of increasing mean v sin i with decreasing 
stellar mass. 



4 DETECTION THRESHOLDS FOR M DWARF STARS 
WITH STARSPOTS 

In this section, we simulate detection thresholds for low-mass plan- 
ets orbiting at the centre of the habitable zone of M dwarf stars. We 
use the starspot models from the preceding section to derive ac- 
tivity induced jitter. Since the parameter space for characterising 
detection thresholds is so large, we have chosen to fix a number 
of parameters as in previous sections. We have simulated the de- 
tection thresholds for 1, 2, 5, 10 & 20 Me planets that orbit in 
the habitable zones of 0.1, 0.2 & 0.5 Mq stars. Table 2 lists the 
periods of the planets for each stellar mass and the corresponding 
stellar radial velocity amplitudes, K„, for each planet. In this sec- 
tion we have simulated all stellar and planetary orbit inclinations 
with i = 90° , and used a range of v sin i values as discussed in 
§3.4. 

For each star/planet combination, we generate radial velocity 
points for a range of observation epochs. We simulate 10, 20, 50, 
100, 200 & 500 epochs, with one observation made every night 
for simplicity. Jitter from the two sources (i.e. starspots and instru- 
mental/measurement precision) is then added to the planetary radial 
velocities. The starspot models 1-6 are used to add the stellar ac- 
tivity jitter to each planetary radial velocity point. This is achieved 
by sampling, at a random observation phase, the line profile that 
arises from the particular starspot model we are interested in. We 
use model estimates to approximate the instrumental/measurement 
precision. This is an important further consideration, since v sin i 
affects the best precision that can be achieved. For a n assumed Y 
band S /N ratio of 90-100 and R~70,000, Fig. 7 of iReiners etall 
(2010) indicates the appropriate accuracy that may be achieved in 
the Y band for a 3000 K atmosphere. This varies from ~ 1.5 m s _1 
at v sin i = 2 km s _1 to ~ 11 m s _1 at v sini = 50 km s~\ An 
example of a simulated planet induced stellar RV signal is shown 
in Fig.|4]before (solid/red) and after (dashed/green) jitter is added. 
This planet is detected with FAP < 0.01 in Fig.|7](upper right panel, 
green square at 100 epochs). 

We then car ry out a Lomb Scargle periodogram analysis 
dPress et alj 1 19921) on each radial velocity curve in an attempt to 
recover the planetary signature. Although this is adequate in our 
simulations, we note that we have not included the effects of or- 
bital eccentricity (i.e. e = in all simulations), which will modify 
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Figure 5. Detection false alarm probabilities (FAPs) vs number of observation epochs for habitable zone planets orbiting M = 0.1 Mq (squares connected by 
solid/red lines) 0.2 Mq (circles connected by dashed/green lines) and 0.5 Mq (triangles connected by dotted/blue lines) stars. The five curves for each stellar 
mass represent planetary masses of 20, 10, 5, 2 & 1 decreasing in a left-to-right sense (labelled for the 0.5 Mq curves). The horizontal line indicates the 
1 per cent false alarm probability (FAP = 0.01) with the grey region representing undetected planets. All points with FAP < 0.01 are considered as detections 
of the planet. The plots are for activity model 2 with v sin i = 2 km s — 1 , instrumental precision = 1.5 m s _1 , and with starspot contrast ratios of T s = 0.65T P 
(high contrast) on the left and T s = T p — 200 K (low contrast) on the right. 
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Figure 6. As for Fig. 5 for activity models 2 (top) and 4 (bottom) with v sin i = 5 km s 1 (instrumental precision = 2 km s ). 



our detection thresholds ( Ford 2 005r) in extreme cases. The results 
of our periodogram analyses are plotted in Figs.[5]-[E]as detection 
thresholds. We have selected to illustrate a range of scenarios that 
represent varying v sin i and activity levels, following the preceding 
discussion. 

The false alarm probability (FAP) is plotted against the to- 



tal number of observations (see caption of Fig. [5] for full de- 
tails). The left hand plots show the results for the cases where 
T a = T p — 200 K (low contrast) while the right hand plots show 
the cases for which T s /T p = 0.65 (high contrast). As expected, 
our ability to detect a planetary signature decreases with increas- 
ing v sin i and stellar activity level. 
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Detection thresholds - model 6, vsini = 10 kms" 1 (low contrast) 
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Figure 7. As for Fig. 5 for activity models 3 (top) and 6 (bottom) with v sin i = 10 km s 1 (instrumental precision = 3 km s 1 ). 



Fig.[5]represents our low v sini and activity (model 2 - active 
solar analogue) level. Since both models 1 & 2 exhibit significantly 
< 1 m s _1 starspot jitter (see Fig. 3) with v sini = 2 km s -1 , we 
illustrate the results for model 2 only since iReiners et all [2010) 
predict that the highest precision achievable in the Y band with 
vsini = 2 km s -1 is ~1.5 m s _1 . The precision therefore domi- 
nates the noise in this scenario rather than the starspot jitter. As few 
as 20-30 epochs of observations are required to detect ^C2 Mq 
planet orbiting a 0.1 Mq star, while 50 epochs are required to de- 
tect a 1 Mq planet. However for a 0.5 Mq star, 500 epochs are 
required to detect a 1 Mq planet. Obtaining such a large number 
of observations is unlikely to be feasible, at least with limited tele- 
scope allocations and a sufficiently large target sample. There is 
little difference between the low and high contrast scenarios, again 
because all starspot induced jitter is well below the achievable pre- 
cision level. 

Fig.[6]illustrates results for v sini = 5 km s _1 and two differ- 
ent activity models (model 2 & model 4 - with up to 9 per cent 
spot filling). In this instance, for model 2, the changes are slight 
compared with the preceding case where v sini = 2 km s~\ This 
observation illustrates a further important limitation to precision, 
namely the instrumental resolution. A greater degree of spotted- 
ness is required (i.e. model 4 cases in Fig. [6} before a noticeable 
change in detection thresholds begins to take place. Here, while 80 
epochs in the low contrast case will enable 1 Me planets to be 
detected, ~500 epochs are required t o dete ct 1 Mq planets in the 
high contrast scenario. IJenkins et all J2009h report a median v sin i 
for 0.5 Mq and 0.2 Mq stars that most closely match the v sin i = 
5 km s _1 used in Fig. [6] This simulation may therefore be deemed 



to represent the detectability of earth-mass to few-earth-mass plan- 
ets around average 0.5 Mq and 0.2 Mq stars. 

In Fig. [7] we show the expected detectability of planets around 
stars with v sini = 10 km s~ . This rotation velocity is the median 
for later M spectral types, such as the 0.1 Mq star. By 10 km s _1 , 
the less active model 3 scenario still enables 2 Mq planets orbiting 
0. 1 Mq stars to be detected in the low contrast case, whereas 5 Mq 
planets may be detected in the same number of epochs for the high 
contrast case. Only planets with mass ^ 10 Mq may be detected 
with < 100 epochs of observations for the high contrast case. With 
a highly active star (model 6), only one planet (20 Mq) remains 
detectable in the high contrast regime, albeit requiring 500 epochs 
of observations. For model 6, the starspot jitter contribution is so 
large (compared with the instrumental/measurement precision) that 
it dominates the noise contribution. We also noted in §3.3 (see also 
Fig. 3) that the jitter is greatest for the 0.2 Mq model with low 
spot contrast spots. As a result, in this scenario, our model predicts 
that detection of a planet orbiting a 0.5 Mq star becomes easier 
than for a planet orbiting a 0.2 Mq star. 

The global trend seen in Figs. 5-7 continues in Fig. [8] where 
only the most massive planets may be detected for very active stars. 
If the spot filling factors of up to 50 per cent, repor ted at earlier 
spect ral types for stars rotating with ~ 20 km s _1 dO'Neal et al.l 
1998) are to be found in M dwarfs, then a lower planet detection 
limit of order 1 Neptune mass may be expected. By 50 km s~ (not 
shown), only 20 Mq planets may be detected in orbit around a 0.1 
Mq star with low contrast spots (model 6) in less than 100 epochs 
of observations (for model 4, only ^10 Mq planets may similarly 
be detected in ^ 100 epochs). 
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Detection thresholds - model 6, vsini = 20 kms" 1 (low contrast) 
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Figure 8. As for Fig. 5 for activity models 3 (top) and 6 (bottom) with v sin i = 20 km s 1 (instrumental precision = 6 km s 1 ). 



4.1 The effects of activity and observing strategy on planet 
detectability 

We have demonstrated how the relative contrast ratio, F p /F 3 is an 
important factor that determines the detectability of planets, par- 
ticularly when spot activity dominates the jitter. However, it must 
also be realised that the observing strategy simulated in this paper 
will determine the detectability of planets. In Particular, the orbital 
period of the planet may be either somewhat shorter, or longer than 
the span of the observations. Table 2 indicates that the orbital peri- 
ods of the simulated habitable zone planets are ~ 5, 13 & 36 days 
for 0.1, 0.2 & 0.5 M© stars respectively. Hence, detection of a pe- 
riodic signal on timescales shorter than the period is less likely. In 
other words, it should be easier to detect habitable zone planets or- 
biting lower mass stars owing to the shorter periods. Since we only 
simulate a minimum of 10 epochs, this mostly applies only to the 
36 day period experienced by a planet orbiting a 0.5 Mq star. It 
can be seen from Figs. 5-8 that no planets are detected in orbit 
about a 0.5 Mq star until over half the period has been sampled. 
Therefore if trageting such stars for habitable zone planets, a differ- 
ent observing strategy would be needed to minimise the number of 
observation epochs. In reality, a radial velocity survey would wish 
to search for planets with a range of orbital radii and periods, espe- 
cially as shorter period planets are more likely to be detected. 

The effect of activity on the number of epochs required to 
detect a planet can also be assessed by completely removing the 
starspot jitter from all simulations. In our low activity simulation 
(vsini = 2 km s _1 instrumental precision of 1.5 m s" 1 and Model 



2 starspot coverage), removing the starspot activity does not have 
any significant effect on the detection curves plotted in Fig. 5. This 
has already been discussed in the preceding section and is a conse- 
quence of the instrumental precision at this rotation velocity (and 
implied instrumental precision) dominating the jitter. The effect of 
removing spot activity from the simulations presented in Figs. 6 - 
8 however leads to significant changes in the number of epochs 
necessary for detection, as might be expected. In our most extreme 
scenario simulated in Fig. 8 (high contrast), although no planets are 
detected, complete removal of the starspot jitter leaves a 6 m s _1 
precision floor due to the rotation velocity of v sini = 20 km s~ . 
Here, planets of 5-20 may then be detected with 300 - 40 
epochs respectively. 



5 SUMMARY & CONCLUSIONS 

We have used more realistic starspot models than previous stud- 
ies to determine the effect of activity induced jitter in precision 
radial velocity studies of M dwarf stars. We have demonstrated that 
with only several tens of epochs, habitable zone earth-mass planets 
can be detected around low-activity stars. Since the contrast ratio 
between photosphere and spots is uncertain, we opted to simulate 
two extreme cases. In reality, the true c ontrast ratio likely falls be- 
tween the two extreme cases simulated. Rockenfeller et al. (2006a) 
required temperature differences of only a few hundred K in order 
to fit their lightcurves of M5 V & M9V stars. However, their mod- 
elling assumed only a single spot was present. A more uniformly 
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spotted surface, especially if highly spotted could produce a sim- 
ilar lightcurve but would be expected to require a higher contrast 
between photosphere and spot to achieve the same photometric am- 
plitude. 

The factors that determine the starspot induced jitter have im- 
portant consequences for the esti mated detection thre shold limits. 
In agreement with other studies l lReiners et alj [2010) we showed 
(Fig. 1) that once the contrast ratio between spot and photosphere 
is sufficiently high, the improvement in jitter, when moving from 
V-band to Y-band, is less impressive. In this case, the spots are sim- 
ply dark enough that the radial jitter does not show such a strong 
decrease with increasing wavelength. Conversely, the decrease in 
jitter as a function of increasing wavelength is much more pro- 
nounced at lower contrast ratios. At low contrast, the exact nor- 
malised line strengths and relative line strengths of the photosphere 
and spot play an important role in determining the exact relative 
jitter at different wavelengths. This additional factor leads to the ir- 
regularity in relative RV amplitudes as a function of increasing T p , 
as shown in Fig. 2 (see §3.2). Although we do not show the full 
jitter amplitude ratio for other photometric bands, Fig. 2 (top) for 
T p = 3250 K indicates that the relationship between wavelength 
and jitter is not as straightforward as a simple blackbody model 
might suggest. Further estimates of the wavelength dependence of 
jitter m ay be obtained by observation of the results o flReiners et al.l 
d2010h who plot RV amplitude as a function of wavelength for the 
5000 A - 18000 A range (their Fig. 12). At longer wavelengths than 
Y band, equivalent width effects at low photosphere/spot contrast 
apart, there is little further gain in precision. A further important 
consideration is the number of lines available for cross-correlation, 
which may be fewer at longer wavelengths. 

We have assum ed that activity scales with rotation 
(Brow ning et alj l2010h so that less spotted stars are slow rota- 
tors and more spotted stars are fast rotators. However, as dis- 
cussed in i]2.1l Reine rs"& Basril j201Ch find that the relation be- 
tween rotation and activity is weaker in M dwarfs. Knowledge 
of the true spottedness of these stars is clearly vital for an accu- 
rate estimation of the detectability limits for low mass planets or- 
biting the lowest mass M dwarfs. In addition, spot coverage de- 
rived from o ther spectroscopic m ethods using temperature sensitive 
(TiO) lines {cy Neal et alj |l998) indicate that there may be large 
discrepancies with results from Doppler images. For active stars, 
up to ~ 50 per cent spot coverage has been derived using the 
TiO method, while Doppler images of similar stars typically derive 
< 10 per cent coverage. The difference between these methods 
is that the TiO procedure derives an unresolved mean spot cover- 
age while the Doppler imaging results are limited by the amount 
of resolvable information in the broadened rotation profile of the 
star. It seems likely that active stars are therefore more spotted than 
Doppler images show, and that they possibly exhibit smaller spots 
than are derived via this method. In summary, high spectral resolu- 
tion Doppler imaging surveys (in the red-optical or infrared), pos- 
sibly combined with other global spot coverage surveys are needed 
for a more accurate picture of spot patterns on mid-late M dwarfs. 

It is clear that there are many physical factors that will deter- 
mine the true detection threshold for planets orbiting M dwarfs; 
we have only simu l ated p lanets in circular orbits for instance. 
I Kennedy & Kenvonl ([2008) simulated planet scattering and found 
that for low mass stars only, planets with long circularisation times 
on eccentric orbits could form. Although we have not simulated 
the effects of eccentricity in full, we have already noted that in ex- 
treme cases, detection t hresholds will be raised and more observa- 
tions will be required jFordll2005h . We find that for longitude of 



periastron, cj = 0°, v sini = 10 km s _1 (instrumental precision = 
3 m s _1 ), the number of nights required to make a detection in- 
creases by approximately 2 - 2.5 times for eccentricity, e = 0.5 and 
by ~ 5 times for e = 0.9. 

In this paper, we have shown that earth-mass or near earth- 
mass planets in the habitable zones of late-mid M dwarfs can be 
detected when the number of starspots matches those seen on the 
Sun at extremes of activity. Moderate rotation quickly increases the 
number of observations that are required to make a detection such 
that by the time vsini = 10 km s _1 , ^ 100 observations are re- 
quired to detect planets with masses ^ 5 . Variations in v sin i 
and contrast ratio are therefore important factors in particular in 
determining the detection thresholds. Given t he evidence for a gen- 
eral increase in v sin i with decreasing mass dDelfosseetaiI l998; 
iMohantv & Basrill2003l ; lienkins et alj|2009h among M dwarfs, the 
balance of these factors may be important. Obtaining estimates of 
starspot coverage from further modelling and photometric obser- 
vations will provide vital information that will enable a link be- 
tween starspots and more traditional chromospheric activity indi- 
cators dPelfosse et al|l99fc ; IMohantv & Basrill2003l ; iReiners et all 
2009i: lBrowning et alfeoiO) to be made. 
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